3 



The asymptotic behaviour of the 
number of solutions of polynomial 

congruences 

Dirk Segers* 
August 23, 2012 

Abstract 

One mentions in a lot of papers that the poles of Igusa's p-adic zeta 
function determine the asymptotic behavior of the number of solutions of 
polynomial congruences. However, no publication clarifies this connection 
precisely. We try to get rid of this gap. 

1 Introduction 

(1.1) Let / e Z[ polynomial over the integers in n variables. Put 

). We are interested in the number of solutions of f(x) = mod m 
in (Z/mZ) n for an arbitrary positive integer m. The Chinese remainder theorem 
reduces this problem to the case that m is a power of a prime. Fix from now 
on a prime p and let Mi, % e Z> , be the number of solutions of the congruence 
f[x) = mod p l in (Z/p*Z) n . The aim of this paper is to study the asymptotic 
behaviour of the numbers M i; and to relate this behaviour with information about 
the poles of Igusa's p-adic zeta function, which will be defined in (1.3). 

(1.2) Let Z p be the ring of p-adic integers. The behaviour of the Mj is well 
understood if nas no singular point in Z™. Indeed, we can take a k G Z >0 
for which / has no singular point modulo p h because / has no singular point 
in the sequentially compact space Zl\ Using Hensels lemma, one obtains that 
Mi = M 2k -iP [n ~ m ~ 2k+1) for every i > 2k - 1. 
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(1.3) The behaviour of the Mi is more complicated if / has a singular point in 
Z™. At this stage, we introduce Igusa's p-adic zeta function Zf(s) of /. It is 
defined by 

Z f (s)= [ \f(x)\ s \dx\ 

for s G C, Re(s) > 0, where \dx\ denotes the Haar measure on Q™, so normalized 
that Zp has measure 1. Note that Zf(s) only depends on p~ s . We will write 
Zf(t) if we consider Zf(s) as a function in the variable t := p~ s . 

All the Mi describe and are described by Zf(t) through the equivalent relations 

Z tlt)= P lt) -?VLzl and P( t ) = l^£W, 
where the Poincare series P(t) of / is defined by 

oo 
i=0 



(1.4) Igusa proved in Igl that Zf(s) is a rational function of p~ s by calculating 
the integral on an embedded resolution of the singularities of /, which always 
exists by Hironaka's theorem [Hi]. This implies that it extends to a meromorphic 
function Zf(s) on C, which is also called Igusa's p-adic zeta function of /. We 
also obtain from the relations in (1.3) that P(t) is a rational function. 

Igusa determined actually a specific form of the rational function which in- 
volves geometric data of an embedded resolution g of /. He obtained that Zf(t) 
can be written in the form 

A(t) 



Z f (t) 



n i6j (l - pr»iW 



where A(t) G S[t], with S := {z/p l \ z G Z, % G Z> }, where A(t) is not divisible 
by any of the 1 — p~ v H Ni and where the Nj and Uj — 1 are the multiplicities of 
fog and g*dx along an irreducible component Ej of g~ 1 (f~ 1 {0}). It is surprising 
that most irreducibele components of 5 ,_1 (/ _1 {0}) do not induce a factor in the 
denominator. This would be elucidated if the monodromy conjecture (see for 
example |Dej ) is true. 

It follows from (1.3) and Zf(t = 1) = 1 that we can write 

nt) - m 



where B(t) G S[t]. Here, B(t) is also not divisible by any of the 1 —p u H N i . The 
poles of P(t) and Zf(t) are actually the same. 
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(1.5) In this paper, we try to explain the relation between the poles (and their 
order) of P(t), which are the same as those of Zf(t), and the numbers Mj. If 
also the principal parts of the Laurent series of P(t) at all poles are known, we 
will even calculate the numbers Mi (and not only their asymptotic behaviour) 
for i large enough. The principal parts of the Laurent series of Zf(t) and P(t) 
at a certain pole can be calculated from each other, which is also the case for 
the ones of Zf(s) and Zf(t) at corresponding poles. Therefore, it is also possible 
to calculate the numbers Mj for i large enough from the principal parts of the 
Laurent series of Zf(s) at all its poles. This will not be worked out in the paper 
because it leads to formulas which are more complicated and which do not give 
us more insight. 

Reference. An introduction to Igusa's p-adic zeta function which contains more 
details can be found in [Sell Section 1.1], |Ig2| or |De] . 

Acknowledgements. I want to thank Pierrette Cassou-Nogues for pointing my 
attention at this problem. 



2 The asymptotic behaviour 

(2.1) We define an equivalence relation on J. We say that ji ~ iff v hl^h = 
Vj 2 /Nj s . This equivalence relation determines a partition of J into sets Jk, k G V. 
For k G V, we denote the lowest common multiple of the Vj, j G Jk, by and 
the lowest common multiple of the Nj, j G Jk, by Remark that cik/bk = Vj/Nj 
for all j G Jk- Let be the cardinality of Jk- Because 1 — p~ ak t bk is a multiple 
of 1 — p~ u H N j for all j G Jk, we can write 

Pit) - °® 



where C(t) G S[t\. 



Theorem. There exists a unique decomposition of every Mi with i > deg(P(t)) 
of the form 

r 

M i = J29k(i)p rihV , 
k=i 

where the Ik are different rational numbers and where every Qk{i) is a nonzero 
function which is polynomial with rational coefficients on residue classes. The r 
numbers lk~n are the real parts of the poles of Zf(s). If we denote the elements 
ofV by 1, ... , r in such a way that Ik — n = —cik/bk for every k G {1, . . . , r}, we 
have for k G {1, . . . , r} that 

1. the function Qk{i) is polynomial on each residue class modulo bj., 
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2. the maximum of the degrees of these polynomials is equal to rrik — 1 and 



3. these polynomials (and thus also gk(i)) are determined by the principal 
parts of the Laurent series of Zf(s) in the poles with real part —dk/bk- 

Remark. (1) The Ik are rational numbers less than n because the real parts of 
the poles of Zf(s) are negative rational numbers. The author proved in [Se3j 
that the real part of every pole of Zf(s) is larger than or equal to —n/2. This 
implies that Ik > n/2 for every k G {l,...,r}. Moreover, in the case that n = 3 
and / has no singular point in l? v of multiplicity 2, the author proved |Se2] that 
there are no poles with real part less than —1, which implies that Ik > 2 for every 
k G {l,...,r}. 

(2) It follows from the theorem that the asymptotic behaviour of the number 
of solutions is determined by the largest real part of a pole of Zf(s) and by the 
largest order of a pole with maximal real part. 

Proof. Applying decomposition into partial fractions in Q[t], we can write 

C k (t) 



P(t)=C (t) + ^-- 



j)-o- k fb k \m k 
k&V x y ' 

where every Ck{t) G Q[t] and where deg(Cfc(£)) < rrikbk for k G V. Note that 
the term C (t) does not give a contribution to Mj for i > deg(C (t)) and that 
deg(Co(i)) = deg(P(t)) if one of them is non-negative. Now we look at the 
contributions of the other terms. So fix k G V. Note that C&(£) contains exactly 
the information of the principal parts of the Laurent series of P(t) at the poles 
with absolute value p ak l bk . We have 

Ck{t) Ck,m k (t) _|_ Ck,m k -l\t) Cfc.lW 



'1 — p-akt b k) m h (1 — p-akfh}™*: [\ — p-a k fb k ^m k -l \ _ p-a k fb k 

6fc— 1 oo 

= J2Y,9Ue)p- eak t ebk+d 

d=0 e=0 

& fc~ 1 °° +eb k +d 

_da k /b k j r (n-a k /b k )(eb k +dy 



d=0 e=0 



pn{eb k +d) 



where Ck,i{t) G Q[t] with deg(Ck,i{t)) < bk and where the maximum of the degrees 
of the polynomials gk,d{e) is equal to rrik — 1. Actually, if we denote the coefficient 
of t d in C k ,i{t) by C k ,i, d , we get 

(e + m fc -l)! (e + m k -2)\ 

9kAe) ~ Lm e\(m k - 1)! + G ^-M e!(mfe _ 2)! + + G ^.i^- 

□ 
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(2.2) Finally, we give two examples. In the first example, all the coefficients 
of the polynomials Ck(t), k G V, are in S. This is in some sense the easiest 
situation. The second example shows that this is not always the case. There are 
several ways to compute the Poincare series: one can calculate the integral on an 
embedded resolution of singularities of /, one can use the formula for polynomials 
which are non-degenerated over ¥ p with respect to their Newton polyhedron [DHJ 
and one can use the p-adic stationary phase formula |Ig2[ Theorem 10.2.1]. All 
these techniques are also explained in |Selj Section 1.1]. 

Example 1. Let f(x, y) = y 2 — x 3 and let p be an arbitrary prime. Then, 
-t 6 + pH 2 - pH 2 + p 6 



P(t) 



{p 5 -t 6 )(p-t) 

2p~H 5 + 2p~H A + 2p~H 3 + 2p- 2 t 2 + (p + l)p- 2 t + (p+ l)p -1 p- 1 



1 - p-H & 1 - p-H 

We obtain for every e G Z>o that 



M 6e = (p+iy*- 1 -p 6e -\ M 6e+1 = ( P +l)p 7e - P 6e , 

M 6e+2 = 2p 7e+2 -p (ie+1 , M 6e+3 = 2 p 7e+3 -p 6e+2 , 

M 6e+4 = 2p 7e+4 -p 6e+3 and M 6e+5 = 2p 7e+5 - p 6e+4 . 



Example 2. Let f(x, y) = x 3 + y 5 and let p be an arbitrary prime. Then, 

-t 15 + {p- l)t 14 + (p - l)pt 12 + (p - l)p 3 t 9 
+(p - l)p 3 t 8 + (p - l)p 5 t 5 + (p - l)p 6 t 3 + (p - l)p 6 t 2 + P 9 
P(t) = (p 8 -t 15 )(p-t) 

C 2 (t) 

l_ p -8 f 15 1-p-l^' 

where 

p 7 +p- 2 14 p 7 +p 2 - p- 1 13 p 7 +p 2 -p- 1 12 p 7 + p 3 -p 2 - 1 n 
11 1 (p 7 -l)p 8 (p 7 -l)p 8 (p 7 -l)p 7 (p 7 -l)p 7 

p 7 + p 3 — p 2 — 1 10 p 7 + p 3 — p 2 — 1 9 p 7 + p 4 — p 3 — l^_ g p 7 + p 5 — p 4 — 1 7 
(p 7 — l)p 6 (p 7 — l)p 5 (p 7 — l)p 5 (p 7 — l)p 5 

p 7 + p 5 — p 4 — 1 6 p 7 + p 5 — p 4 — 1 5 p 7 + p 6 — p 5 — 1 4 p 7 + p 6 — p 5 — 1 3 
+ (p 7 - l)p 4 t + (p 7 - l)p 3 1 + (p 7 - l)p 3 * + (p 7 - l)p 2 1 
2p 7 - p 6 - 1 2 p 8 - 1 ^ p 8 - 1 
(p 7 — l)p 2 (p 7 — l)p 2 (p 7 — l)p 



' 2(t) " 'W^Wp' 
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As an illustration, we calculate the Mj for i in the residue class of 3 modulo 15: 

M 3 + 15e - (1) 

= p 4 + 22e + (P - l)/ +22e (P"l)p 2+15e 



P 7 — I p 7 — 1 

p 4 +22e+(p _ 1) P + 2+ 15e 

P' — 1 

p 4 +22 e + (p _ 1)(p 7e + . . . + p U + p 7 + 1)p2+ l 5e> 



Note that the two terms in ([T]) are not integers and that one of them is negative. 
Note also that the Poincare series in the two examples are rational functions of t 
and p, but this is not the case in general. 
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